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Introduction 

Entanglement is one of the interesting features of quantum systems. It is used as a physical 
resource in realization of some important quantum information and quantum computation 
processes such as quantum parallelism [T], quantum cryptography [2J, quantum teleportation 
[3l [5] , quantum dense coding [61 [7] , reduction of communication complexity |8j and beating 
classical communication complexity bounds with entanglement [9]. In these applications usu- 
ally a source produces entangled particles and after receiving these particles to related parties, 
there is an important question for the parties - are these particles already entangled? 

One approach to distinguish entangled states is entanglement witness (EW) [9l [10] • A 
quantum state is entangled iff there exists a Hermitian operator W with Tr(Wp) < and 
Tr(Wp sep ) ^ for any separable state p sep [13]. In this case, the witness, W, detects the 
entanglement of the density matrix, p. Entanglement witnesses have some advantages with 
respect to other efficiently implementable witnesses such as concurrence [T2]. Recently there 
has been an increased interest in nonlinear EWs because of their improved detection with 
respect to linear EWs. A nonlinear EW is any bound on nonlinear function of observables 
which is satisfied by separable states but violated by some entangled states [9]. 

Optimization problems occur in both classical and quantum physics [19]. One of the im- 
portant subclasses of these problems, is convex optimization where the related functions of the 
problems are convex. The importance of convex optimization is that in these optimizations, 
any locally optimal solution of the problem is guaranteed to be globally optimal [2D] . On the 
other hand, the set of all possible states of a quantum system that can occur in nature must 
be a convex set [2 1 j and as the state of a quantum system is fully characterized by density 
matrix of that system so the density matrix must be a convex function. Therefore, convex 
optimization is a natural optimization in quantum information. 

In this article we discuss the problems and advantages of the convex optimization approach 
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for finding entanglement witnesses of a given system. The approach considered here, is similar 
to that of article |14j . with some improvements such as feasible region for one particle, which 
this one enable us for detecting entanglement of bipartite systems in d dimension. In particular, 
it will be shown that for some density matrices the whole entanglement region can be detected. 
This approach is essentially goal-directed but it may not work for all possible density matrices. 

Organized into five sections, the first section reviews the definition of an entanglement 
witness for a quantum system. The next section begins with the feasible region definition for 
a single qudit followed by some feasible regions example for qudits. A technique of convex 
optimization algorithm for finding EWs is presented in section 3. Section 4 focuses on some 
examples for 3x3 bipartite systems. Finally, section 5 deals with a class of d x d density 
matrices, including examples for entanglement detection conditions for d — 3, d — 4, and 
d = 5. 

1 Entanglement Witness 

A bipartite mixed state is called separable, if it can be prepared by two parties in a classical 
way, that is, by agreeing over the phone on the local preparation of states. A density matrix 
that has been created in this way can only contain classical correlations. Mathematically this 
means: a mixed state p is called separable iff it can be written as p2] 



otherwise it is entangled. Here the coefficients pi are probabilities, i.e. ^ Pi ^ 1 and 



An entanglement witness W is a Hermitian operator such that Tr(Wp s ) ^ 0, for all sepa- 
rable states p s , and there exists at least one entangled state p e which can be detected by W, 
that is Tr(Wp e ) < 0. The existence of an EW for any entangled state is a direct consequence 
of Hahn-Banach theorem [T5] and the fact that the space of separable density operators is 
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convex and closed. 

For a bipartite d x d system, one can expand an EW as follows 

1=1 i=i 

where Id is a d x d identity matrix, A^j are some parameters, and O^s ( Ojs ) are iVi ( A^ 2 ) 
Hermitian operators from first ( second ) party Hilbert space. 



2 Single particle FR 

Consider the following mapping 

p s i-> Tr(Wp s ) (2.3) 

This is a mapping from convex region of separable density matrices to a convex region called 
feasible region (FR). As a separable density matrix, p s , can be written as a convex combination 
of product states ( 11. ip . for one party ( say first one ) this map has terms like 

\a i ){a i \^Tr{6 i \a i ){a i \) = P i (2.4) 

The set of PjS defines a convex region for first particle and we will refer to this region as single 
particle FR. 

In this section we determine single particle FR. To this aim consider a single particle, a 
d-level quantum mechanical system, or a qudit. The possible states of this particle can be 
represented by state vectors in the associated d-dimensional Hilbert space. Each observable 
associated with this particle, is represented by Hermitian or self-adjoint linear operators acting 
on the corresponding Hilbert space. These operators can be expanded in terms of the traceless 
Hermitian matrices which are the matrix representations of d 2 — 1 independent infinitesimal 
generators of SU(d) group ( the special unitary group which is a matrix Lie group of dimension 
d 2 — 1). So in the following we use these generators as EWs operator expansion basis. Specially, 
we consider two sets of operator expansions for a d level single particle as follows. 
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In this set we choose d 2 — 1 operators 



where 7j's, the infinitesimal generators of SU(d), are represented as traceless Hermitian matri- 
ces for a d dimensional system. The FR for this set of operators is defined by ( see Appendix 
A for a proof ) 

d 2 -i 

(2-6) 

Second set 

In this set we choose d 2 operators as follows. Every d-dimensional square matrix can be 
written in terms of square matrices E.., which has the value 1 at the position and zeros 
elsewhere. Now one can define Hermitian traceless basis for d-dimensional matrices as follows 
[T6] . The off-diagonal basis are given by 

O a ,f>=-j=(E atf , + E Pia ), a<(3 
O a ,fs = -^=(E af3 -E 0a ), a>(3 

and the diagonal basis are 



Oa,a E aa 

where a, (3 = 1, ...,d and E a ^ = \a)(/3\. For this set, the FR becomes ( Appendix A ) 

d 2 

a,/3=l 

3 Convex optimization formalism for bipartite systems 

In this section we will express the construction of an EW for an arbitrary bipartite state in 
terms of a convex optimization. A convex optimization problem [17] , consists of optimizing an 
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objective under some linear or nonlinear matrix equality and inequality constraints, precisely, 
we deal with a problem as follows 

minimize fo{x) 

subject to fa(x)^0, i = l,...,m. (3-8) 
hj(x)=0, j = l,...,p. 
Here, the objective and the constraint functions are convex functions which means that they 
satisfy inequality fa (ax +(3y) ^ afa(x)-\-(3fa(y), for all x, y G R and alia, (3 G .R with a+/3 = 1, 
a ^ 0, (3 ^ and the equality constraint functions hi(x) = must be affine (A set C G R n is 
affine if the line through any two distinct points in C lies in C). 

We present a two-step convex optimization algorithm as follows. In the first step, recalling 
the definition of an EW, we impose the first condition, Tr(Wp s ) ^ 0, for all separable states 
p s . As a matter of fact, in this step Tr(Wp s ) is the objective function and inside of the FR 
will be defined by the inequality constraints. In the second step, for a given p, we impose 
the second condition for an EW, Tr(Wp) < 0. Now, the objective function ( which will be 
minimized ) is Tr(Wp), and the inequality constraints come from the solution of the first step. 
The following sub-sections provide more details on these two steps. 

3.1 First step, the nonnegativity conditions for Tr(Wp s ) 

Here the approach is similar to [18] but for the convex optimization. Consider a non-positive 
Hermitian operator of the form in (11.21) . We will attempt to choose the real parameters 
Aij such that W becomes an EW. To this aim, we introduce the maps Pi = Tr(Oip s ) and 
Pj = Tr(0'jp s ) for any separable state p s . These maps, map the convex set of separable states 
into a bounded convex region which will be called feasible region. The first property of an EW 
is that its expectation value over any separable state is non-negative, i.e., the condition 

Tr(W Ps ) = l + J2A, j P i P->0 
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must be satisfied for any point (Pi,P 2 , ...) of the feasible region. In order to satisfy this 
condition, it is sufficient that the minimum value of Tr(Wp s ) be non-negative. Using standard 
convex optimization, we find this minimum value and then we impose the non-negativity 
condition on this minimum. 

For the bipartite system, as we studied before, the FRs are defined by ££i Pf < 1 and 
Ylfjli Pf ^ 1- These constraints in matrix notation become P*P < 1 and P ft P' ^ 1, where 
the vector P is P = (P i; P 2 _ i )*. Then the problem can be written as a convex optimization 
problem in the form 

Minimize 1 + P l AP' 
subject to P*P ^ 1 

, P' l P' <: l. 

The Lagrangian associated with this problem, is given by 

L = 1 + P'AP' - -\{P l P - 1) - -A'(P'*P' - 1) (3.9) 

Any points satisfying K arush-Kuhn- Tucker (KKT) conditions, are primal and dual optimal, 
with zero duality gap. These conditions are 

1. primal constraints: P*P - 1 and P n P' - 1 < 0, 

2. dual constraints: A >z and A' ^ 0, 

3. complementary slackness: A(P*P — 1) = and A'(P' P' — 1) = 0, 

4. gradient of Lagrangian with respect to P, P' vanishes: VL = 0. 

Requiring these conditions for our problem, the forth condition VL = reduces to 

AP' = AP (3.10) 

and 

P t A = \'P' t (3.11) 
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and the complementary slackness condition requires that A(F*P-1) = and \'{P' P'-l) = 0, 
which reduce to P t P = 1 and P' P' = 1 noting that our problem is a minimization problem, so 
A > and A' > 0. Multiplying (13.101) by P l from left side and using the reduced complementary 
slackness condition, P l P = 1, yields 

P f AP' = \P l P = A 
also multiplication (13. lip by P' from right side and using P ' P' = 1, leads to 

P l AP' = X'P'tp' = A' 
then A = A'. So, from (13. lOj) and (13. lip one can write 

A l AP' = \ 2 P' 

and 

AA l P = X 2 P 

where, \f for i — 1, d 2 — 1, are eigenvalues of A 1 A. If we choose A's in a way that Af ^ 1 
(i.e. —1 ^ | Aj| ^ 1) for all z's, then the minimum of Lagrangian (13. 9p . 1 + Al, is nonnegative 
which leads to the nonnegativity of Tr(Wp a ). 

3.2 Second step, the detection condition for a given density amtrix 

For a given density matrix p g , the entanglement detection condition is the negativity of 
Tr(Wp g ) so we want to minimize this term. Again this problem can be written as a con- 
vex optimization problem in the following form 

Minimize Tr(Wp g ) 

s.t. A 1 A -I d ^0 
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As we will see, the solution to this problem yields a nonlinear EW. Lagrangian, L, associated 
with this problem is 

L = Tr(Wp g )+Tr\(A t A-I d )Z] 

where Z is a positive matrix which plays the role of the Lagrange multipliers associated with 
the inequality constraint A 1 A — I d < 0. 

Requiring the KKT conditions for this problem, as in the previous subsection, yields the 
equations 

p + A(Z t + Z) = 

and 

A 1 A = I d , 

where p is a matrix with components 

p^^Tripfi^O^). (3.12) 
If Z be a symmetric matrix (i.e., Z t = Z) then p = —2AZ so that p l p = 4Z 2 or 

Z = ±(p t p)h (3.13) 

and 

A = -\pZ-\ (3.14) 

The matrix A, in general, has a nonlinear form which yields a nonlinear form for the associated 
EW. Now the minimum of the Lagrangian becomes 

l-Tr[y/pfl (3.15) 

and the negativity of this term is the entanglement detection condition for the given density 
matrix. 

The remaining part of the paper gives examples for two-qutrit systems and higher dimen- 
sional systems. As we will see in these examples, this method mainly gives necessary and 
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sufficient condition(s) for separability. However, this method can not detect entanglement for 
all bipartite systems. In [22], the authors have provided a general method to improve the en- 
tanglement detection for bipartite systems via nonlinear expressions in EWs. Their approach 
is an iterative method and can yield useful approximation. However, there is only one example 
for two-qubits, and no example for higher dimensions. 



4 Examples for qutrit systems 



A. The two qutrit Horodecki states 

As first example illustrating the methodology, consider the state described in [23] given by 



2 a 5 — a 



(4.16) 



with < a < 5, |V>+> = tiELN), ^+ = |(|01><01| + |12)(12| + |20)(20|) and a_ = 
|(|10)(10| + |21)(21| + |02)(02|). This state is PPT for 1 < a ^ 4 and as it was shown in 



the state is bound entangled for 3 < a ^ 4, separable for 2 ^ a ^ 3 and free entangled 
for 4 < a ^ 5. In [24], Doherty et al. introduced a non-decomposable EW which detects the 
Horodecki bound entangled states for 3 < a ^ 4. 

Here we work with the second set of operators introduced in the second section, the EW 
takes the form 

9 

W = h ® h + Ai >A ® °i 
and the Z matrix in (13.131) is equal with 



126 















7 + 2c 
7 - c 
7 - c 










7 - c 
7 + 2c 
7 - c 










7 - c 
7 - c 
7 + 2c 
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where, c = a/19 + 3a (a — 5). The detection condition is given by 

1 - Vl9-15a + 3a 2 < 

which detects all of the entanglement region, i.e. for ^ a < 2 and 3 < a ^ 5. So the 
necessary and sufficient condition for separability is obtained. 

B. Density matrix in unextendible product bases (UPBs) 

The density matrix considered here, is constructed using UPBs [25] which is a bound 
entangled density matrix for 3 <8> 3 state. The state has the following matrix form 



/ 



1 

72 



V 



-7 
-7 

2 

2 

2 

2 

2 

2 

2 



-7 
-7 

2 

2 

2 

2 

2 

2 

2 



2 
2 

-7 

2 
2 

-7 
2 
2 
2 



2 
2 
2 
2 
-16 
2 
2 
2 
2 



2 

2 

2 

2 

2 

2 

2 
-7 
-7 



2 
2 
2 
2 
2 
2 
2 

-7 
-7 



(4.17) 



Using the second set of FR relations, the detection is Tr(Wp L 



A is 



V 



C. Horodecki state 



0.2909 
0.1587 
0.2909 







-0.8127 
0.3316 
0.1873 



0.1587 
0.8463 
0.1587 







0.3403 
-0.0416 
0.3403 



0.2909 
0.1587 
0.2909 







0.1873 
0.336 
-0.8127 



0.1873 
0.3403 
-0.8127 







-0.2868 
-0.1551 
-0.2868 



0.3316 
-0.0416 
0.3316 







-0.1551 
-0.8545 
-0.1551 



—0.087 and the matrix 
\ 



-0.8127 
0.3403 
0.1873 







-0.2868 
-0.1551 
-0.2868 



(4.18) 



In third example, we consider the 3x3 bound entangled state introduced in Ref. [26] as 



Pa 















l + q 
2 





(4.19) 
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This state is PPT while it is entangled for all < a < 1. Working with the second set of 
operators, we plot the detection in terms of parameter a, which shows that, all PPT region 
has been detected by our method ( see Figure 1 ). 



5 Generalized qudit Choi maps 



In order to improve our knowledge of the convex optimization formalism for detecting entan- 
glement of quantum systems, we will concentrate here on exploring d x d bipartite Choi maps 
density matrices [28]. These density matrices have the following form 

1 - p d_1 1 
Pppt =/#o,o)(V'o,o| + Yl Pipu ( 5,2 °) 

i=l 

where |Vo,o) = ^ E?=o 1010, < p < 1, Pi = E?=o 10(^1 ® \l + *|, » = 0, 1, d - 1 

and fiiS are parameters with properties ^ ^ 1, Y^t=x fa = 1- Also ^ nas snown that ( 
the so called PPT conditions ) 



p TA >0= > p< mm{mm{ _^ ; _JWk + / }}, i^zj^zk— l,...,d— 1. 



1 + jUi ' /ij/ife - 1 Y (jUjjUfc - l) 2 ' 

(5.21) 

In the following we study the entanglement detection problem by the convex optimization 
formalism discussed previously and we give detection conditions in general. 

For density matrices (15.201) . working with the second set of operators, we have 

hs,^ = Tr(p PPT O af3 ® Op,), a,P,n,u = 0,...,d-l. 
The off diagonal elements of this matrix are 



Pao.« = — J^Ve-*' a^/3 = 0,...,d-l, (5.22) 



where (3 — a is calculated in modulo d. The diagonal elements are 



Po-.aa = |. a = 0,...,d-l l (5.23) 
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also, 



and finally for a > (3 we have 



P M = j> a< P' 



P 



p 



The relations (15.221) and (15.231) form the following d x d block sub-matrix 

d-i 



= ti 1 ~ p 

Pdxd d d d 



where S is the shift operator for d ® d defined as 



S 



V 



i=l 

as 



o) 



1 
1 



The eigenvalues of (p 4 p) d xd are 



1 



-b+(l-p)T4|, fc = 0,...,d-l 



where 14 = J2j=i Pj^ k an d w = e Finally the detection condition (13.151) for total density 
matrix, (15.201) . is equal with 

d-l 



k=0 



< 0. 



(5.24) 



Now in the following we study the cases d = 3 and d = 4 in details, and we will obtain the 
necessary and sufficient conditions for separability in d = 3. 

A. The case d=3 



It is instructive to examine our formalism for the case d = 3 which leads to necessary and 
sufficient conditions for the entanglement detection. For d = 3 explicitly we have 



1 -p 1 

Pppt = P|V'o,o)\'0o,o| H g— (A*iPi + A*iPi), P^g' Pi + P2 = 1 



(5.25) 
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where |^,o) = ^(|0)|0) + |1)|1) + |2)|2)) , Pl = |0)(0| ® |1)(1| + |1)(1| ® \2)(2\ + |2)(2| <g> |0)(0| 
and p 2 = 1 0> <0 1 ® |1)(1| + |1)(1| ® |2)(2| + |2)(2| <g> |0)(0|. The PPT condition for this density 
matrix is given by 



^ p < 



A*iA*2 



Using the second set of operators, the detection condition for this PPT density matrix (I5.24p 



is 



2 - 6p - [(1 - 3p) 2 + 3(p - 1) 2 (1 - 2/n) 2 ] 2 < 0. 
Both PPT and detection conditions are satisfied for the following regions 



2 1 + n\ 



/il /Xi^2 

< P ^ 7 + 



/Xl/i 2 



P>\P>i 



(/xi/i 2 - l) 2 



or 



2 2-/ii /ii/i 2 - 1 Vl^i^-l) 



Other PPT regions, ^ p ^ •j+ 3 7 an d ^ V ^ iT^T' wn i cn could not be detected, are cor- 
respond to separable regions ( for a proof see Appendix B ). Therefore applying our convex 
optimization formalism for density matrices (15.251) yields the necessary and sufficient condi- 
tions for separability 



B. The case d=4 

For d = 4 the entanglement detection condition (15.241) reduces to 



1 3 



k=0 



< 0. 



which is equal with 



1 - [3p + -(1 + 2\p + (1 - p)V!\ + \p + (1 - p)V 2 \ 



(5.26) 



where we use V = P i + /i 2 + /a 3 = 1 and V\ = V 3 . Simplifying this relation and using 
Hi = 1 — /ii — /i3 yields 

^{3 - 12p - \p + (1 -p){l - 2m - 2fi a )\ 
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- 2vV - 2p{l - p)(l - /n - /i 3 ) + (1 - P ) 2 [l - 2/^(1 - /n) - 2/4,(1 - /is)]} < 0. (5.27) 
We can remove the absolute value signs and replace that terms with the following ones 

p + (l-p)(l-2 f i 1 -2fi 3 ) (5.28) 

and 

-\p + (1-^(1-2//! -2fi 3 )) (5.29) 

So, the detection condition (I5.27P with PPT condition (15.211) . define two possible detection 
regions. These regions are presented in tablel and table2 of the appendix C, which are corre- 
spond to (15.281) and (I5.29P respectively. In Figure 2. these regions are plotted. 



C. The case d = 5 

For d — 5 the detection condition (I5.24p is equal with 

1 - {Ap + ^[1 + 2(\p + (1 - p)Vx\ + \p + (1 - p)V 2 |)]} < 0. (5.30) 

As a special case we assume ImV\ = which gives fxi = |[(1 — Vo)/^ + (Vo — l)//3 + 2/i 4 ] 
and /i2 = ^ 3 ( 1+v/ ^j^ 2 ^ 4 ~ 1 ) where we use Y^t=i fa = 1- The detected PPT region under these 
simplifying assumptions plus PPT conditions (15.211) for d — 5, is plotted in Figure 3. 



Summary and conclusion 

In this paper we have defined single particle feasible region ( FR ) and have used it for 
constructing entanglement witnesses for mixed bipartite systems by the convex optimization 
method. The advantage of considering the single particle FR is that it produces a straightfor- 
ward formulation in the convex optimization method. We have also given an explicit method to 
construct nonlinear EWs for given bipartite systems. The ability of this approach is illustrated 
by applying it to the separability problem for some well known two qutrit systems. However, 
this method can not detect entanglement for every bipartite systems. 
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This approach leaves some questions open such as the optimality of these EWs. So a natural 
next step would be to investigate the optimality of the obtained EWs. Is there any relationship 
between these EWs ( specially, those EWs leading to necessary and sufficient conditions for 
separability ) and the optimal entanglement witnesses condition considered in [3U] ? However, 
such problems which are under investigation, perhaps will probe some aspects of optimal EWs 
problems. 
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Appendix A 

a). Proving (12.6ft 

First we prove that working in any basis, will not change the quantity X^Li 1 p f- Recalling 
Pi mapping definition, Pi = Tr{O i \a){a\) 1 we have 

d 2 -i d 2 -i 

i=l i=l 

If we change the basis |a) = U\f3) where U is a d x d unitary operator, then 

d 2 -i <f-i 

J2 p " = J2^ u1 °^^ u1 °^ 

i=l i=l 

Using the adjoint representation of the Lie algebra, 

d 2 -i 

3=1 

where Oj's are generators of the Lie group, and A's are a subgroup oidxd orthogonal matrices. 
In fact one can expand the generators of the Lie algebra as above. So we have 

d 2 -l d 2 -l d 2 -l d 2 -l 

E^ 2 = E(E A ^m°m^\°^) = E m°^)\ 2 

i=l j,k=l i=l i=l 

( the term in parenthesis, J^f-^ 1 is equal to 5j^) which completes the proof. 

Now we choose \/3) = [1,0,0, ...,0] T , a d — 1 dimensional vector and a appropriate set of 
generators, 7,, [27J, 

&j,k = Pj,k ~t~ Pk,j 

Vj,k = -i(Pj,k ~ p k,j) 
1 



Wi 



1(1 



k=l 

where P m ,n = 1 ^ j < k ^ d and 1 ^ I ^ d — 1. All of these operators are traceless 

and they generate the algebra su(d). With these choices, we obtain 

E^ 2 = Ek/w>i 2 = (^n) E i^i^i^i 2 

i=l i=l ^ > i=l 
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All of the operators in ( lA-il) vanish in summation unless w t s, which for them we have 



and as the number of such terms is d — 1 so 



:^ry)ElW'')l 2 = (^T))E I (T T T; 



d ,44,2 2 s d 2, 



*2(d-l) y |^Z Z + l^ ( 2(d-l) )(2 d) 1 



As we interest to the entire FR, then 



d 2 -i 



i=l 
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Appendix B 

a). Proving separability 

Here we prove that for PPT regions ^ p ^ and ^ p ^ <^T> ^ ne density matrix 
(I5.25P is separable. 

For first region, ^ p ^ i+t^j if p = then p PPT = ^(pipi + pipi) which is obviously 
separable and if p = then p PPT could be written as p PPT = y^-{3/Ui[|(|^o,o)(V'o,o| + ^ + 
y)] + (1 — 2yU!)y}. The term [|(|^o,o)(V'o,o| + ^ + y)] has been explicitly represented as a 
mixture of product states [29J so the later p PPT is separable. 

For second region the separability proof is similar to the first region. 
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Detected entangled regions for generalized qudit Choi maps, d — 4 
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Table 1. 



A*3 


/il 


P 


di < fJ, 3 < ^ 




Ci < p ^c 2 


b 2 < fii < b 3 


C3 < p < c 2 


Hi = &3 


P = c 3 


4 < A*3 ^ 7 


bi < fii < b 2 


Cl < p ^ C2 


b 2 < /ii < & 3 


C3 < P < c 2 


i < ^3 < | 


6i < Ui < /!•? 


Ci < » ^ Co 


U-? < Ui < 69 


Ci < 10 ^ C9 


Wi = 69 


< 10 ^ C9 


b 2 < fii < b 3 


C3 ^ P ^ c 2 


Ail = &3 


P = c 3 


= \ 

t*>6 3 


64 < ui < 4 


C4 < » s> c^ 


| < A*3 < § 


61 < Ail < &2 


Ci < p < c 2 


/il = &2 


C3 < P < c 2 


6 2 < Ail < & 3 


c 3 < P < c 2 


Ail = ^3 


P = c 3 




Aii = h 


c 3 < P < c 2 


Aii = h 


P = c 3 



where 

ai = ^(9 - 4^2), a 2 = 2^(9 + 4v^), 

61 = 4 + 97/i 3 - 28^3 + 12a4 h = £ I 48 " 73 ^3 - 4^(3^3 - I) 2 ], 
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6 3 = |(4 - 7/x 3 + ^3(8-15^)), h = |(109 - 28^/15), 



_ 28 + 6(^ +/i 3 ) / 15(/xg + (4) - 12(/n + /i 3 ) + 6/X1/X3 + 4 

Cl (/i 1 -/i 3 ) 2 -16(/i 1 + /i 3 )-32 V [(//! - // 3 ) 2 - 16(//! + // 3 ) - 32] 2 

/ ^i/j 3 ~ _|_ 1^3 _ 2(/ii+^ 3 )-l _ 3^i(3^i-6y / l5-32)+6v / 15-65 _ ygi 

2 V (M1M3-1) 2 M1M3-1' 3 2(^1+^3) ' ° 4 3/xi(3/ii-50)-335 ' ° 5 v^ZT+v^' 

Table 2. 



A*3 


A*i 


P 


< /i 3 < a 3 


ni < /ii ^ 65 


c 6 < p < c 2 


65 < nx < 6 6 


c 6 < p < c 7 


^3 ^ A*3 < a l 


n 2 < /ii < 65 


c 6 < p < c 2 


&5 < /"I < &6 


c 6 < p < c 7 


ai < /x 3 < 1/3 


&5 < A*l < ^6 


c 6 < p < c 3 


&2 < Aii < 63 


c 6 < p < c 3 


63 < /ii < 67 


c 6 < p < c 2 


1/3 < /i 3 < a 2 


^2 < A*l < ^3 


c 6 < p < c 3 


63 < /ii < 65 


c 6 < p < c 2 


&5 < A*i < ^8 


c 6 < p < c 7 


A*3 = a 2 


ai < a*i < &9 


C 8 < p < Cg 


69 < /ii < 610 


C 8 < p ^ C10 


a 2 < /i3 < 1 


n 3 < /ii < 65 


c 6 < p < c 2 


65 < /ii < 6 8 


c 6 < p < c 7 



where 

_ 2641-1740^ 



6 5 = i[2 - ^3 - ^3(4-3^)], 6 6 = 1 - 2^3, &r = |(2 - At 3 ), &s = |(1 - A* 3 ), 
6 9 = ^(49 - 4^2 - a/705 + 248^), &10 = |(5 - >/2), 
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. , 1 + 6f4 - 3^3 + 6/i§ - 3/ii(l + /i 3 ) 22 - 3/ii - 3//. 

C 6 = 1 -2W — — -t ; ; - ; ; r + 



7 M1+M3-2' 



[/i 2 - 2 Atl ( At3 - 4) + ^ 3 (8 + - 32] 2 /1 2 - 2^/13 - 4) + /i 3 (8 + ^ 3 ) - 32 

( 29(12^2- 611 + 87 fix) 

° 8 ~ + 24711 - 1000^ - 29^(214 - 8^ + 29^ 

-58* 



736 + 84V2 + 174/xi(29/ii - 19 - 2y/2) 



[lOOOv/2 - 24711 + lOOOv/2 + 29/ii(214 - 8^ + 29^i)] 2 



r - 1 1 / 29(9+4y / 2) m 29 _ i , 29_ 



[(9+4>/2)/ii -29] 2 9+4^)^1-29' iu 4^-49+29^1' 

Although iii,ni,ni exist but they have not simple algebraic form so we did't write them. 
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Figure Captions 

FIG. 1: The entanglement detection in terms of a. All entangled states are detected for 
a =[0,1]. 

FIG. 2: The entanglement detection region for the d = 4 case in terms of parameters /ii 
and /i 3 . All states in this region are PPT entangled states. 



FIG. 3: The entanglement detection region for the d = 5 case in terms of parameters /i 3 
and (14. All states in this region are PPT entangled states. 
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